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SECTION  1 


INTRODUCTION 

Radio-wave  scintillation  has  been  used  extensively  to  study  irregu¬ 
larities  in  a  variety  of  transmission  media;  moreover,  with  continued 
Improvements  in  technology  and  data  analysis  procedures,  additional 
details  of  the  irregularity  characteristics  are  coming  into  focus.  Yeh 
and  Liu  [1982]  reviewed  the  Interpretation  of  ionospheric  scintillation 
data,  which  is  our  Interest  here.  The  essential  features  of  the  theory 
can  be  Illustrated  by  considering  the  simplest  model — a  thin  scattering 
layer.  Under  weak,  scattering  conditions,  the  average  spectrum  of  the 
intensity  scintillation  measured  in  a  plane  perpendicular  to  the  propa¬ 
gation  diectlon  admits  the  representation 

♦  j(K)  =  4r^X^AL*jj  (K,0)sln^(K^t^)  ,  (1) 

”e 

where  r  is  the  classical  electron  radius,  X  is  the  wavelength,  AL  is 
e 

the  layer  thickness,  and  ♦  (K,K  )  is  the  three-dimensional  spectrum  of 

N  Z 

e 

the  electron-density  structure  evaluated  with  K  =  0  where  K  lies  along 

z  z 

the  propagation  direction,  z.  The  Fresnel  radius,  l^,  is  defined  as 


=  /xz/(4i) 


(2) 


References  are  listed  in  alphabetical  order  in  Section  4. 


2  2  2 

The  term  sin  (K  i^)  In  Eq.  (1)  substantially  alters  the  otherwise 

direct  mapping  of  the  In  situ  structure  as  characterized  by  ♦  (K,K  ), 

N  Z 

e 

particularly  for  spatial  frequencies  smaller  than  2i/t^.  If  the 
structure  Is  uniformly  distributed  along  an  extended  path  but  still 
produces  only  weak,  scattering,  there  Is  an  additional  filter  function 
that  acts  to  limit  the  depth  of  the  K-dependent  oscillations.  Thus,  the 
pronounced  diffraction  nulls  predicted  by  Eq.  (1)  are  rarely  observed  In 
actual  data  F Crane.  1976],  but  the  high-pass  filter  characteristic  Is 
preserved. 

For  phase  scintillation,  the  representation  corresponding  to  Eq. 

(1)  is 

♦.(K)  =  rfx^AL*„  (K,0)  cos^(K^tJ)  ,  (3) 

♦  ®  "e  ^ 

which  Is  similar,  but  there  Is  no  attenuation  for  K  <  l/l^.  Thus,  the 
envelope  of  «^(K)  Is,  potentially,  a  more  direct  mapping  of  the  In  situ 

9 

Irregularity  structure  than  the  corresponding  envelope  of  ♦^(K).  None¬ 
theless,  the  Intensity  spectrum  generally  provides  the  most  consistent 
measure  of  scintillation  strength  through  the  scintillation  Index  S^, 
defined  by  the  relation 

=  <I^>/<I>^  -1  =  II  ♦^(K)dK/(2T)^  .  (4) 

—00 

The  Intrinsic  suppression  of  (K,0)  In  ♦_(K)  for  K  <<  separates 

N  I  f 

e 

trend-llke  variations  from  the  more  homogeneous  random  structure.  Thus, 
as  long  as  Is  less  than  unity.  It  can  be  used  as  a  measure  of  the 
overall  perturbation  level  [Rlno  and  Liu.  1982]. 

We  can  only  measure  a  time  series  generated  by  the  combined  motion 

of  the  source  and  the  Irregularities.  The  envelope  of  ♦.(K)  cannot  be 

9 

obtained  directly:  however.  If  the  Irregularities  do  not  change  over  the 
time  of  the  measurement,  a  t ime-to-sp>ace  translation  can  be  specified  In 


The  definition  of  T  depends  on  n  FRino.  1982],  but  Eq.  (8)  shows  that 
the  power-law  index  of  the  Intensity  power  spectriun  beyond  the  Fresnel 
cutoff  is  simply  related  to  the  one-dimensional  spectral  index,  n. 
Similarly,  the  phase  power  spectrum  in  a  single  power-law  environment 
satisfies  the  relation 


(9) 


I 

I 

I 


•'ll 

i 


I 


I 

i:? 

I* 


<P,(f)  a  Tf 


essentially  over  the  entire  range  of  f;  moreover,  the  fundamental 
features  of  these  results  are  preserved  In  an  extended  medium — although 
the  relation  between  T  and  becomes  more  complicated. 

In  previous  analyses  of  scintillation  data  from  the  Wideband  satel¬ 
lite  FFremouw  et  al..  1978;  Rlno,  1979;  Fremouw.  1980;  Livingston  et 


al. .  1981;  Secan  and  Fremouw.  1983; 


1985],  the  model 


outlined  above  was  used.  A  large  body  of  data  was  summarized  by  per¬ 
forming  log-linear  least-squares  fits  to  the  measured  phase  spectra  over 
fixed  temporal  frequency  ranges.  As  described  In  Livingston  et  al. 
[1981]  and  Rlno  [1979],  both  the  equatorial  and  auroral-zone  data  showed 
fairly  shallow  and.  In  the  case  of  the  equatorial  data,  systematically 
varying  spectral  Indices.  The  main  characteristics  of  the  equatorial 
Irregularities  have  been  confirmed  by  In  situ  measurements,  which  also 
revealed  a  more  complicated  two-component  form  for  the  In  situ  irregu¬ 
larity  structure  [Livingston  et  al. .  1981;  Rlno  et  al. .  1981a;  Basu  and 


Basu.  1983; 


:,  1986].  More  recent  analyses  of  other 


scintillation  data  have  also  shown  evidence  of  two-component  spectra 
[FranRe  and  Liu.  1983;  Bhattacharwa  and  Rastogl.  1985].  Comparable 
analyses  of  auroral-zone  data  have  Just  begun,  but  It  appears  that 
similar  results  apply  [Villain  et  al . .  1986]. 

In  light  of  these  new  results,  it  was  appropriate  to  reanalyze  at 
least  a  representative  subset  of  the  equatorial  and  auroral-zone  data 
from  the  Wideband  satellite  experiment.  We  looked  for  evidence  of 
multicomponent  spectra  and  compared  the  data  with  those  that  are  cur¬ 
rently  being  acquired  with  the  HILAT  satellite  (which  was  launched  into 
an  82  degree  Inclination  circular  orbit  at  830-lun  altitude  In  June 
1983).  VHF  and  UHF  amplitude  and  phase  scintillation  data  from  the 
HILAT  satellite  have  been  collected  at  ground  stations  In  Sondre 
Stromfjord  (Greenland),  Tromso  (Norway),  Fort  Churchill  (Canada),  and 
Bellevue  (Washington),  but  we  shall  consider  only  the  data  from  the 
auroral-zone  station  at  Tromso  in  this  comparative  study. 


Ideally,  to  perform  such  an  analysis,  we  would  use  numerical 
computations  with  an  hypothesized  spectral  form  in  the  weak-scatter 
theory  and  would  compare  the  predictions  to  the  measured  spectra.  This 
requires  extensive  computation,  however,  and  it  is  practical  only  for 
small  amounts  of  data.  Where  multifrequency  data  are  available,  as  with 


Wideband  and  HILAT,  the  measured  power  spectra  can  be  tested  for  adher¬ 


ence  to  the  frequency  dependence  predicted  by  the  weali-scatter  theory. 
For  example,  Eq.  (8)  is  frequency  independent,  and  Eqs.  (7)  and  (9) 


scale  Inversely  with  frequency  squared.  The  Fresnel  nulls  themselves, 
which  depend  on  wavelength  and  distance  through  Eq.  (2),  can  also  be 
measured.  In  a  companion  report  FLlvlngston  and  Dabbs.  1986  (hereafter 


referred  to  as  Report  I)],  this  procedure  was  used  both  to  establish  the 


limits  of  the  weah-scatter  theory  and  to  determine  the  contaminating 


effects  of  the  diffraction  nulls  In  measuring  the  characteristics  of 


multicomponent  spectra,  it  was  found  that  the  most  accurate  results  are 


obtained  by  fitting  Individual  power-law  segments  adjacent  to  the  first 


diffraction  null. 


It  Is  desirable,  however,  to  have  an  automated  procedure  that  can 


be  routinely  applied  to  the  very  large  data  bases  that  have  already  been 


and.  In  the  case  of  HILAT.  continue  to  be  accumulated.  To  this  end,  a 


curve-fitting  procedure  has  been  developed  that  will  search  for  a 


specified  maximum  number  of  power-law  segments.  Both  the  power-law 


segments  and  the  break  frequencies  are  determined  automatically  and 


consistently.  Because  the  curve-fitting  procedure  has  a  variety  of 


applications,  it  is  described  In  detail  in  the  Appendix.  In  the 
remainder  of  this  report,  we  describe  (1)  the  results  of  applying  the 


automated  curve-fitting  procedure  to  the  old  Wideband  equatorial  and 


auroral -zone  data,  and,  (2)  In  the  latter  case,  the  comparisons  to  more 


recently  acquired  HILAT  auroral-zone  scintillation  data.  The  indepen¬ 


dent  analysis  of  the  equatoral  data  summarized  In  Report  I  provides  a 


basis  for  Identifying  any  contamination  caused  by  departures  from  strict 


power-law  behavior,  which  is  induced  by  the  diffraction  nulls. 


SECTION  2 


DATA  ANALYSIS 


2.1  GENERAL  APPROACH. 

From  the  discussion  In  the  Introduction,  it  is  clear  that  even  in 
the  case  of  weak,  scattering,  it  is  difficult  to  fully  accommodate  all  of 
the  effects  that  can  influence  the  Interpretation  of  the  scintillation 
spectra  in  terms  of  the  structure  of  the  in  situ  irregularities.  On  the 
other  hand,  certain  characteristics  of  the  measured  spectra  seem  to  have 
a  robust  relation  to  the  corresponding  in  situ  characteristic.  For 
example,  the  power-law-index  relation  given  by  Eq.  (8)  does  not  depend 
on  the  geometry  or  *^.  Of  course,  if  the  local  value  of  the  spectral 
index  (or  any  other  characteristics  of  the  in  situ  structure)  varies 
along  the  propagation  path,  we  will  infer  a  weighted  average  value,  but 
this  problem  is  Inherent  in  any  integral  measure,  and  it  must  ultimately 
be  resolved  by  detailed  modeling  and  by  comparisons  with  in  situ 
measurements,  which  are  being  pursued  separately. 

In  our  analysis  here,  we  accept  the  fundamental  limitation  of  the 
Integral  measure  and  look  for  evidence  of  multicomponent  power-law 
structures.  Intuitively,  we  expect  that  as  long  as  each  power-law  seg¬ 
ment  subtends  a  large  enough  range  of  scale  sizes,  relations  such  as 
Eqs.  (7)  and  (9)  can  be  applied  to  each  power-law  segment.  This  is 
discussed  in  Rlno  and  Owen  [1982],  and  has  been  demonstrated  in  numer¬ 
ical  simulations  (Franke  and  Liu.  1  83;  Rino  and  Owen.  1984].  Analytic 
relations  between  parameters  that  characterize  the  three-dimensional, 
multicomponent  in  situ  spectra  and  the  measurable  spectral-strength 
parameters  are  more  illusive,  but  good  approximations  are  beginning  to 
emerge  FWittwer .  1986]. 
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In  our  data  analysis,  we  use  the  intensity  scintillation  level  to 
ensure  that  we  are  within  the  confines  of  the  weak-scatter  theory  as 
summarized  In  Section  1,  and  we  loolc  for  evidence  of  multicomponent 
spectra  In  the  phase  scintillation  data.  Using  the  Wideband  satellite 
data,  phase  scintillation  was  measured  simultaneously  at  VHF  (138  MHz), 
UHF  (413  MHz),  and  L  band  (1239  MHz).  At  UHF,  seven  closely-spaced  fre¬ 
quencies  were  available.  These  data  have  been  used  to  establish  the 
maximum  scintillation  level  for  which  the  weak-scatter  theory  applies, 
as  described  In  Report  I.  The  data  are  processed  In  overlapping  30-s 
segments;  both  linear  detrending  and  windowing  are  applied  prior  to 
standard  spectral  analysis  using  digital  Fourier  transforms.  The 
spectral  estimates  are  then  Interpolated  onto  a  uniform  logarithmic 
grid,  and  the  curve-fitting  procedure  described  In  the  Appendix  Is 
applied.  Thus,  our  summary  data  consist  of  a  number  of  power-law 
segments  that  are  characterized  Individually  by  a  strength  measure,  a 
power-law  Index,  and  a  break  frequency. 


The  limitations  of  the  curve-fitting  procedure  Itself  are  discussed 
in  the  Appendix.  Ultimately,  however,  we  rely  on  the  Internal  consis¬ 
tency  of  the  results  against  our  understanding  of  the  propagation  theory 
and  their  relation  to  other  results  to  verify  the  numerous  manipulations 
performed  In  the  course  of  analyzing  the  data.  Our  principal  concerns 
in  this  regard  are  the  potential  biases  induced  by  diffraction  effects 
and  the  departures  from  the  weak-scatter  theory — both  of  which  are 
discussed  In  Report  I.  We  shall  refer  to  any  functional  dependence  on 
as  a  diffraction  effect. 

2.2  EQUATORIAL  DATA. 

The  general  morphology  of  scintillation  occurrence  as  deduced  from 
Wideband  satellite  data  has  been  summarized  by  Livingston  [1980].  For 
our  analysis  here,  19  passes  from  the  Kwajaleln  station  and  24  passes 
from  the  Ancon,  Peru,  station  were  selected.  This  same  data  set  was 
used  In  previous  studies  [Rlno  and  Owen.  1981,  1982;  Rlno  et  al. . 

1981b].  The  Individual  passes  were  drawn  from  active  nighttime  periods 


but  encompass  a  wide  range  of  activity  levels  and  elevation  angles.  The 
passes  are  Identified  In  Table  1.  Figure  1  summarizes  the  distributions 
of  Intensity  and  phase  scintillation.  The  ordinates  measure  the  percent 
occurrence  of  the  abscissae,  which  are  Identified  on  the  Individual 
histograms.  In  Figure  2  we  show  scatter  plots  of  the  measured  UHF  and 
VHF  Indices  plotted  against  the  logarithm  of  the  L-band  Index  (top 

row  of  plots)  together  with  scatter  plots  of  the  measured  rms  phase  at 
adjacent  frequencies  (bottom  row  of  plots). 

Several  points  should  be  noted.  From  Figure  1,  we  see  that  the 

scintillation  levels  reach  saturation  (S.  ~  1.0)  at  VHF.  moderate  values 

4 

at  UHF,  but  are  generally  weak,  at  L  band.  Thus,  most  of  the  VHF  data 
cannot  be  analyzed  using  the  weak-scatter  theory.  The  saturation  of  the 
VHF  values  Is  more  clearly  displayed  In  Figure  2;  moreover,  a  careful 
Inspection  of  the  data  shows  that  a  small  subset  of  the  VHF  S^  values 
exceed  unity.  This  so-called  strong-focusing  effect  In  power-law 
environments  Is  Indicative  of  an  outer-scale  cutoff  or  contributions 
from  steeply-sloped  power  laws  (Rlno  and  Owen.  1984;  Jakeman .  1982]. 
Evidently,  It  was  not  observed  In  our  earlier  analyses  because  of  the 
shorter  data  Intervals  (10  s  versus  30  s)  (Rlno  and  Liu.  1982].  Turning 
to  the  phase  data,  a  necessary  condition  for  negligible  diffraction 
effects,  as  we  defined  them  above.  Is  a  linear  dependence  of  the  rms 
phase  on  frequency.  The  L-band  and  UHF  data  follow  the  linear  relation 
to  nearly  5  rad  of  rms  phase,  but  the  corresponding  VHF-UHF  data  In 
Figure  2  show  significant  departures,  which  are  indicative  of  strong- 
scatter  effects  that  are  not  accommodated  In  the  theory  outlined  in 
Section  1. 

In  our  summary  analysis  of  the  phase  spectra,  we  consider  only  data 
segments  for  which  the  value  at  the  particular  frequency  lies  between 
0.2  and  0.5.  The  lower  limit  Is  used  to  avoid  noise  contamination  In 


Table  1.  Summary  of  equatorial  passes 


File 

- 1 

Day 

Start 

(0) 

End  (UT) 

Number 

of  Records 

KWAJALEIH 

1977 

(Midnight  = 

12:00  UT) 

188 

12: 

13 

32 

12 

25 

02 

47 

189 

12: 

53 

34 

13 

04 

04 

43 

209 

12: 

13 

28 

12 

23 

58 

43 

212 

12: 

27 

12 

12 

38 

42 

47 

5 

216 

13: 

23 

12 

13 

31 

17 

44 

6 

218 

12: 

57 

12 

13 

07 

42 

43 

224 

11: 

42 

32 

11 

53 

17 

44 

226 

11: 

18 

42 

11 

27 

12 

35 

226 

13: 

01 

52 

13 

12 

07 

42 

10 

228 

12: 

36 

22 

12 

47 

52 

46 

11 

229 

11: 

32 

42 

11 

42 

42 

41 

12 

229 

13: 

17 

02 

13 

26 

17 

37 

13 

232 

11: 

46 

52 

11 

58 

07 

46 

14 

234 

11: 

22 

52 

11 

32 

07 

38 

15 

235 

12: 

01 

22 

12 

12 

52 

47 

16 

237 

13: 

21 

52 

13 

30 

22 

35 

17 

238 

12: 

15 

52 

12 

27 

07 

46 

18 

241 

12: 

30 

52 

12 

42 

22 

47 

19 

250 

11: 

31 

29 

11 

41 

59 

43 

AHCOH  1977  (Midnight 

*  4 

00  UT) 

1 

40 

4 

40 

52 

4 

52 

22 

47 

2 

42 

4 

15 

52 

4 

27 

22 

47 

3 

49 

3 

41 

32 

3 

53 

02 

47 

4 

63 

4 

15 

12 

4 

26 

42 

47 

5 

74 

4 

35 

22 

4 

46 

52 

47 

6 

75 

3 

31 

52 

3 

41 

52 

41 

7 

76 

4 

10 

22 

4 

21 

52 

47 

8 

78 

3 

46 

02 

3 

57 

02 

45 

9 

79 

4 

25 

12 

4 

36 

42 

47 

10 

80 

3 

22 

12 

3 

31 

12 

37 

11 

80 

5 

05 

32 

5 

15 

47 

42 

12 

81 

4 

00 

32 

4 

12 

02 

47 

13 

82 

4 

40 

02 

4 

51 

32 

47 

14 

83 

3 

36 

22 

3 

46 

37 

42 

15 

83 

5 

21 

02 

5 

30 

02 

37 

16 

85 

4 

55 

12 

5 

06 

12 

45 

17 

88 

3 

26 

42 

3 

36 

12 

39 

18 

88 

5 

10 

22 

5 

20 

22 

41 

19 

89 

4 

05 

02 

4 

16 

32 

47 

20 

90- 

4 

44 

42 

4 

56 

12 

47 

21 

96 

3 

31 

02 

3 

41 

17 

42 

22 

103 

4 

39 

12 

4 

50 

42 

47 

23 

104 

3 

35 

32 

3 

46 

17 

44 

24 

110 

4 

04 

22 

4 

15 

52 

47 

Frequency  comparisons  of  intensity  and  phase  scintillation. 


the  midrange  of  the  temporal  frequencies  that  contribute  to  the  power 
spectra.  The  upper  limit  was  established  by  performing  a  detailed 
analysis  of  the  closely-spaced  UHF  scintillation  data,  as  discussed  in 
Report  I. 

The  scintillation  at  each  frequency  is  characterized  separately. 

To  convert  the  temporal  frequencies  to  spatial  frequencies,  we  use  the 

effective  scan  velocity  as  defined  in  Rino  [1982]  for  an  axial  ratio  of 

100:1.  The  basic  assumption  here  is  that  the  anisotropy  of  the  in  situ 

spectrum  can  be  characterized  by  a  quadratic  form — irrespective  of  the 

overall  spectral  shape.  For  equatorial  data,  the  strict  validity  of 

this  assumption  is  not  too  critical.  Indeed,  near  the  geomagnetic 

equator,  Is  nearly  equal  to  the  component  of  the  relative 

penetration-point  velocity  perpendicular  to  the  propagation  path  for 

axial  ratios  greater  than  10:1.  In  Figure  3,  we  have  plotted  the 

distributions  of  v  for  the  equatorial  passes.  The  relatively  small 
GZ  z 

values  of  compared  with  the  actual  penetration-point  velocity  of 

~  3  lun/s  result  from  the  highly  Inclined,  sun-synchronous  orbit  of  the 
Wideband  satellite,  which  places  the  descending  passes  nearly  along  the 
geomagnetic  meridian. 

To  characterize  the  weak-scatter  phase  spectra  in  terms  of  the 

power-law  segments  found  by  applying  the  curve-fitting  procedure 

described  in  the  Appendix,  we  use  the  occurrence  of  breaks  in  the 

middle-frequency  range  (0.32  Hz  <  f  <  3.2  Hz).  Short  power-law  segments 

that  lie  outside  this  range  are  ignored.  Using  this  procedure,  all  of 

the  measured  spectra  were  found  to  be  one  of  three  types:  (1)  single- 

component  power  laws,  (2)  two-component  power  laws  with  p^  <  p^,  and  (3) 

two-component  power  laws  with  p  >  ?„•  In  Figures  4,  5,  and  6,  we 

A  D 

summarize  the  spectral  characteristics  of  the  subset  of  the  phase 

scintillation  data  segments  at  each  frequency,  that  satisfy  the  weak- 

scatter  criteria  0.2  <  <  0.5. 
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The  first  columns  in  Figure  4  contain  histograms  of  the  values 
that  satisfy  the  weak-scatter  limits.  The  VHF  data  have  the  smallest 
populations  because  most  of  the  data  are  too  strongly  disturbed; 


moreover,  the  L-band  data  are  typically  near  the  lower  limit.  The 
center  columns  contain  the  distributions  of  the  spectral  types  as 
described  above,  and  the  right-hand  columns  contain  the  distributions  of 
the  spectral  index  values  for  the  single  power-law  spectra  that  we  have 
referred  to  as  Type  1. 

Using  this  sorting  criteria,  we  find  that  more  than  half  of  the 
data  segments  admit  a  single  power-law  characterization.  The  average 
power-law  index  value  is  near  but  somewhat  greater  than  3.  The 
Kwajaleln  data  give  slightly  shallower  p  indices  than  the  Ancon  data, 
but  the  Kwajaleln  data  have  a  higher  noise  floor  because  a  smaller 
antenna  was  used.  Thus,  we  suspect  that  noise  is  biasing  the  Kwajaleln 
data,  particularly  at  VHP,  and  the  differences  are  probably  not  signifi¬ 
cant.  The  fact  that  the  single  power-law  segments  are  not  confined  to 
the  weaker  absolute  perturbations  is,  however,  significant.  Indeed,  the 
fraction  of  the  samples  that  admits  a  single-component  characterization 
is  essentially  the  same  at  each  frequency.  In  the  in  situ  rocket  data 
summarized  in  Section  1,  single-component  power-law  spectra  were  ob¬ 
served  primarily  in  regions  of  low  electron  density,  such  as  within 
depletions  or  at  higher  altitudes.  Unless  such  structures  are  dis¬ 
tributed  along  a  highly  extended  propagation  path,  we  would  not  expect 
them  to  be  associated  with  the  strongest  scintillation. 

The  characteristics  of  the  two-comp>onent  Type  2  spectra  are 

summarized  in  Figure  5.  The  first  columns  show  the  distributions  of  p^. 

The  center  columns  show  the  distributions  of  break  frequencies  converted 

to  length  units  using  v  The  right-hand  columns  show  the  dlstrlbu- 

er  I 

tions  of  p  .  The  corresponding  data  for  the  Type  3  spectra,  which  are 
B 

included  primarily  for  completeness,  are  shown  in  Figure  6.  Indeed,  the 
number  of  samples  of  Type  3  spectra  is  sufficiently  small  that  we  sus- 
p)ect  their  occurrence  to  be  spurious.  Thus,  we  shall  only  summarize  the 
detailed  structure  of  the  Type  2  spectra. 

We  first  note  that  p.  is  generally  less  than  3,  while  p_  is  greater 

A  B 

than  3;  moreover,  the  distributions  of  p^  show  less  spread  than  the  p^ 

distributions.  As  with  the  p  distributions,  the  p  distributions  from 

B 


the  Kwajalein  data  are  slightly  shallower  than  the  corresponding  distri¬ 
butions  from  the  Ancon  data.  We  attributed  this  to  noise  contamination, 
and  it  follows  that  the  higher  sensitivity  Ancon  data  are  most  accu¬ 
rately  reflecting  the  high-frequency  and  single  power-law  spectral-index 
values;  however,  the  details  of  the  data  selection  and  the  curve-fitting 
procedures  have  subtle  effects  on  the  results.  The  sorting  criteria 
used  here  combined  with  the  tendency  of  the  curve  fitter  to  respond  to 
abrupt  changes  outside  the  middle-frequency  range  evidently  maximize  the 
number  of  Type  1  spectra.  The  procedure  used  in  Report  I,  for  example, 
produced  somewhat  larger  percentages  of  Type  2  spectra. 

The  distributions  of  the  break,  points  shown  in  Figure  5  encompass 
the  values  reported  from  in  situ  rocket  F Rino  et  al . .  1981a;  LaBelle  and 
K ;v.  1986]  and  satellite  FBasu  and  Basu.  1983]  measurements;  however, 
there  is  a  systematic  Increase  in  the  break  scale  in  going  from  UHF  to  L 
band.  If  this  were  truly  reflecting  a  structure  change,  it  would  be  in¬ 
dicative  of  a  perturbation  strength  dependence  of  the  break  frequency, 
as  discussed  below.  Report  I,  however,  shows  that  the  first  Fresnel 
null  in  the  phase  scintillation  data  falls  very  near  the  break  fre¬ 
quency,  and  is  undoubtedly  biasing  the  measurement.  Thus,  we  believe 
that  the  variation  of  the  break  frequency  is  a  bias  such  that  the  L-band 
value  is  closest  to  the  in  situ  value.  The  measured  spectral  character¬ 
istics  are  summarized  in  Table  2  together  with  the  number  of  samples  for 
each  spectral  type. 

Because  we  have  based  our  data  selection  on  the  S,  levels  at  each 
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frequency,  we  have  effectively  sorted  the  data  by  perturbation  strength. 
The  L-band  data  that  exceed  the  lower  limit  arise  from  the  largest 
path-integrated  in  situ  perturbation  levels;  the  UHF  data  are  next.  The 
VHF  data  that  fall  below  the  upper  limit  correspond  to  the  smallest 
path-integrated  in  situ  perturbation  levels.  The  single-component 
spectra  has  almost  no  variation  with  the  frequency  at  which  the  spectra 
are  measured.  This  can  be  seen  more  clearly  in  Figure  7  where  we  have 
plotted  the  measured  VHF  and  L-band  p  indices  against  the  corresponding 


Table  2.  Summary  of  spectral  characteristics  derjved  using 
multicomp)onent  curve-fitting  procedure. 


Frequency 

Type 

No.  of  Data 
Samples/Total 

Spectral 

Index 

(P  or  P^,Pgl 

Break. 

(m) 

KWAJALEIN 

L  Band 
VHF 
UHF 

L  Band 
VHF 
UHF 

L  Band 


27/44 

— 

3.02 

194/295 

3.24 

72/109 

3.16 

9/44 

2.29,  2.96 

66/295 

2.23,  3.18 

25/109 

2.04,  3.43 

8/44 

— 

35/295 

— 

12/109 

— 

ANCON 

74/134 

3.18 

190/345 

3.46 

111/185 

3.36 

39/134 

2.21,  3.35 

142/345 

2.24,  3.54 

70/185 

1.86,  3.43 

21/134 

— 

13/345 

— 

4/185 

— 

These  results  are  subject  to  small  biases  caused  by  the  residual 
diffraction  null  that  falls  near  the  break  frequency.  Based  on  the 
analysis  in  Report  I,  we  believe  that  the  most  representative  median 
values  for  the  Type  2  spectra  are  p^=  2.3,  p_=  3.9,  and  X  =  650  m. 


(a)  KWAJALEIN 


L-band  S,  value,  which  as  we  discussed  In  Section  1,  Is  proportional  to 
4 

the  path-integrated  perturbation  level. 


The  Type  2  p.  Indices,  however,  show  a  distinct  decrease  In  going 
A 

from  UHF  to  L  band.  In  Figure  8,  the  Individual  data  points  are  shown 


together  with  the  log-linear  least-squares  fits. 
[1981],  using  In  situ  satellite  data,  and  Rlno  e\ 


[1981a],  using  In 


situ  rocket  data,  observed  a  systematic  variation  of  the  In  situ 
spectral  Index  that  would  manifest  Itself  In  the  phase  spectrum  as 


=  Pq  -  nlOlogj^pT  ,  (10) 

where  T  Is  the  power-law  strength  parameter.  Their  value  of  n  a  0.03  Is 

close  to  the  average  values  for  the  two  fits  shown  In  Figure  8.  Because 

the  break,  frequencies  themselves  vary  systematically  with  frequency, 

however.  It  Is  clear  that  the  residual  of  the  diffraction  nulls  Is 

biasing  the  breakpoint  estimates  and  the  spectral  slopes  as  well.  Thus, 

as  discussed  In  Report  I,  If  we  adhere  strictly  to  the  weak-scatter 

criteria,  the  possibility  that  the  p.  trend  In  the  phase  data  Is  a 

A 

diffraction  effect,  cannot  be  excluded. 

There  Is,  however,  other  evidence  of  a  p^  dependence  on 

A 

perturbation  strength.  From  this  same  data  set,  Rlno  and  Owen  [1982, 
198^]  showed  that  by  using  Eq,  (10)  the  measured  coherence  time  for  the 
strong-scatter  VHF  data  could  be  brought  Into  excellent  agreement  with 
predictions  based  on  the  strong-scatter  theory.  Similarly,  In  Rlno  et 
al.  [1981b],  It  was  shown  that  the  behavior  of  the  frequency  decorre¬ 
lation  under  strong-scatter  conditions  Is  also  consistent  with  a  varying 
spectral  Index.  Unfortunately,  the  strength  measures  for  these  analyses 
were  derived  from  the  phase  scintillation  data  rather  than  from  the  L- 

band  S,  Index,  which  Is  a  more  robust  measure  of  overall  perturbation 
4 

strength  [Rlno  and  Liu.  1982].  Moreover,  as  discussed  In  Section  1, 
recent  analyses  of  hlgh-resolutlon  In  situ  data  from  the  AE-E  satellite 
do  not  show  the  systematic  variation  that  was  found  In  the  low- 
resolution  data  and  In  the  PLUMEX  rocket  data  [ Livingston .  private 


communication,  Rlno  et  al. .  1981a].  This  suggests  that  the  trend  Is 
confined  to  spatial  wavelengths  larger  than  5  kjn. 

For  reference,  in  Figure  9,  we  have  plotted  and  p  versus 

using  the  thin  phase-screen  theory  for  a  single-component  power  law,  and 

the  p  variation  implied  by  Eq.  (10)  with  an  assumed  100-k.m  layer 

thlcfcness.  The  plot  verifies  the  nearly  linear  relation  between  log 

and  log  for  a  single-component  power  law  having  a  systematically 

varying  spec‘-ral  index — which  is  necessary  if  we  are  to  use  at  L  band 

as  a  measure  of  absolute  Integrated  perturbation  strength.  For  a  two- 

component  power  law.  a  frequency-dependent  correction  factor  must  be 

added  that  acts  to  reduce  the  value  of  C  .  That  is,  a  given  value  of  C 

s  s 

will  produce  a  smaller  value  than  would  be  inferred  from  Figure  9; 
however,  in  light  of  the  most  recent  analyses  of  high  resolution  AE-E 
data,  this  phenomenon  clearly  needs  further  study. 

2.3  AURORAL-ZONE  DATA. 

The  general  morphology  of  auroral-zone  scintillation  as  measured  by 
the  Wideband  satellite  is  desct ibed  in  Rino  and  Matthews  [1980].  Some 
unique  features  attributed  to  the  distribution  of  structure  and  the 
anisotropy  of  the  auroral-zone  irregularities  are  described  in  Fremouw 
[1977];  Rino  et  al.  .  [1978];  Rino  and  Owen.  [1980];  and 
Livingston  et  al . .  [1982].  For  comparison,  we  have  selected  25  HILAT 
passes  recorded  at  Tromso,  Norway,  which  has  nearly  the  same  Invariant 
latitude  as  the  Wideband  station  at  PoRer  Flat,  Alaslta.  The  HILAT 
satellite  carries  a  modified  version  of  the  Wideband  beacon  with  trans¬ 
missions  at  the  same  VHF  and  UHF  frequencies;  however,  the  L-band  signal 
on  HILAT  's  used  for  telemetry  and  as  a  phase  reference.  Thus,  no  L- 
band  scintillation  data  are  available.  The  Individual  passes  are 
summarized  in  Table  3. 
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Table  3.  Summary  of  auroral-zone  passes. 


Day  I  Start 
POKER  FLAT 

59  13:15 

60  12:10 

61  f:21 

61  11:05 

61  12:50 

62  8:17 

62  1C:01 

62  11:45 

62  13:29 

63  8:56 

63  10:40 

63  12:25 

64  9:36 

64  11:20 

68  8:46 

68  10:30 

82  11:03 

82  12:48 

86  10:13 

86  11:58 

88  9:48 

88  11:33 

88  13:17 

99  10:07 

99  11:52 

102  12:06 


End  (UT) 
(Midnight 


3b 

00:43:27 

00:53:57 

36 

C0:17:23 

00:27:53 

37 

01:33:20 

01:43:50 

38 

01:07:18 

01:17:48 

39 

00:41:16 

00:51:46 

40 

00:15:11 

00:25:41 

40 

23:49:07 

23:59:37 

43 

00:39:06 

00:49:36 

44 

00:13:01 

00:23:46 

44 

23:46:57 

23:57:42 

46 

22:54:55 

23:05:25 

48 

00:10:49 

00:21:19 

48 

23:44:50 

23:55:20 

49 

23:18:45 

23:29:30 

50 

22:52:43 

23:03:28 

51 

22:26:41 

22:37:11 

52 

23:42:38 

23:53:08 

54 

22:50:32 

23:01:17 

56 

21:58:29 

22:08:59 

57 

,  i: 32:25 

21:43:55 

58 

22:48:22 

22:58:52 

59 

22:22:21 

22:32:51 

60 

21:56:17 

22:07:02 

61 

21:30:15 

21:40:45 

62 

21:04:08 

21:14:38 

Number  of  Records 
10:00  UT) _ 

42 
47 
47 
47 

44 
38 
47 
47 
41 

45 
47 

46 

47 
47 

44 
47 
47 

45 
47 
47 
47 
47 

43 
47 
47 
47 


TROMSO  1984  (Midnight  =  : 00  UT) 
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The  display  format  Is  the  same  as  for  the  equator  data.  In  Figure 

10  we  show  histograms  of  the  occurrence  of  S,  and  o^,  and  in  Figure  11 

4  4^ 

we  show  scatter  diagrams  of  the  individual  measurements.  The  scintilla¬ 
tion  levels  are  much  lower  than  the  equator  data,  so  that  the  UHF  and 
VHF  data  comprise  most  of  the  data  that  fall  within  the  weak-scatter 
window.  Also,  the  Tromso  station  uses  a  simple  helix  antenna,  whereas  a 
small  UHF/L-band  dish  and  a  four-element  VHF  helix  array  were  used  with 
the  earlier  Wideband  data.  Thus,  the  Tromso  data  admit  more  noise  and 
multipath  contamination.  Finally,  the  Tromso  data  were  recorded  near 
solar  minimum,  while  the  earlier  Wideband  data  were  recorded  during  a 
period  of  increasing  solar  activity.  In  spite  of  these  differences, 
however,  the  general  occurrence  patterns  at  the  two  auroral-zone 
stations  are  very  similar. 

Because  of  the  complicated  anisotropy  of  the  nighttime  auroral-zone 

irregularities  rLlvlngston  et  al..  1982],  the  conversion  from  temporal 

to  spatial  frequency  is  not  as  straightforward  as  it  is  near  the 

geomagnetic  equator,  where  the  diffraction  pattern  is  essentially  one- 

dimensional.  Thus,  in  Figure  12,  we  show  the  distributions  of  v  for 

ef  f 

10:1  rod-like  irregularities  and  10:5:1  sheet  or  wing-like  irregulari¬ 
ties.  The  orientation  of  the  second  irregularity  axis  is  along  the 
local  L  shell.  The  differences  are  significant.  We  have  therefore 
presented  our  summary  results  using  two  estimates.  We  note  that 

the  Wideband  satellite  was  in  a  sun-synchronous  orbit,  whereas  the  HILAT 
satellite  orbit  drifts  through  24  hours  of  local  time  in  6  months.  The 
geometry  of  the  HILAT  passes  relative  to  the  local  magnetic  meridian 
will  vary  more  than  the  corresponding  geometry  of  the  Wideband  satellite 

passes.  Nonetheless,  the  v  distributions  for  10:5:1  wing-like 

ef  f 

irregularities  are  very  similar  for  the  two  satellites. 

Using  the  same  spectral  characterization  procedure  as  that  used  for 
the  Equatorial  data,  we  obtain  the  distributions  shown  in  Figure  13.  No 
L-band  data  from  the  Wideband  satellite  survive  the  low  limit,  which 


TYPES 


is  consistent  with  the  generally  lower  perturbation  levels  that  are 
characteristic  of  the  auroral  zone.  Unllk.e  the  equator,  in  the  auroral 


zone,  Type  3  spectra  with  p^  <  p^  constitute  the  dominant  two-component 
configuration.  The  single-component  p  indices  are  slightly  shallower 
than  the  corresponding  population  of  Equatorial  passes,  but  are  other¬ 
wise  similar.  The  lower  sensitivity  of  the  Tromso  data  clearly  mani¬ 
fests  itself  in  the  distribution  of  p,  compared  with  the  corresponding 
Wideband  data,  but  the  general  characteristics  of  the  spectra  obtained 
at  the  two  stations  are  remarkably  similar. 

Because  of  the  very  small  number  of  Type  2  spectra  in  the  auroral 

zone,  the  distributions  are  not  well  defined,  and  we  have  not  considered 

these  data  further.  The  Type  3  spectra  are  summarized  in  Figure  14.  We 

see  that  the  p.  distributions  are  virtually  identical  between  the  two 
A 

auroral-zone  stations,  having  a  median  value  between  3  and  4.  The 
steeply-sloped  portion  of  the  spectrum  extends  at  least  to  one 
kilometer.  The  variation  is  due  to  the  different  anisotropy  models  used 
to  infer  the  actual  break  point,  as  discussed  above.  Measured  aniso¬ 
tropy  parameters  are  not  available  for  these  data.  The  high-frequency 

component  of  the  Type  3  spectra  is  very  shallow  with  p  approaching  2; 

B 

however,  the  high-frequency  ends  of  the  spectra,  even  in  the  Poker  Flat 

data,  are  susceptible  to  noise  contamination  at  the  lower  S,  levels — 

4 

these  results  are  likely  to  be  biased  toward  lower  values.  The  general 
characteristics  are  consistent  at  the  two  stations.  The  median  values 
are  summarized  in  Table  4. 

The  shallowly-sloped  p  component  of  the  Type  3  auro:'al-zone 

B 

spectra  is  similar  to  the  shallowly-sloped  p^  component  found  in  the 

equator  data  at  somewhat  larger  scale  sizes  and  possibly  caused  by  a 

similar  process;  however,  within  the  limits  of  the  available  data,  we 

cannot  establish  a  perturbation  strength  dependence  of  p  in  the 

B 

auroral-zone.  In  the  earlier  analysis  of  the  Wideband  satellite  data. 


(a)  L  BAND 


Table  4.  Summary  of  auroral-zone  spectral  characteristics 
derived  using  multicomponent  curve  fitter. 


No.  Of  Data 

Spectral 

Break,  (m) 

Frequency 

Type 

Samoles/Total 

Index 

10:1:1 

10:5:1 

POKER  FLAT 

VHF 

r 

1 

214/508 

3.10 

UHF 

1 

74/159 

3.12 

— 

VHF 

2 

26/508 

UHF 

2 

11/159 

— 

VHF 

3 

268/508 

3.44,  2.11 

1072 

912 

UHF 

3 

74/159  j  3.46,  2.19 

813 

398 

TROMSO 

VHF 

1 

258/529 

1 

2.99 

— 

UHF 

1 

87/237 

2.60 

— 

VHF 

2 

201/529 

— 

— 

UHF 

2 

15/237 

-- 

— 

VHF 

3 

251/529 

3.39,  1.18 

660 

309 

UHF 

3 

135/237 

3.23,  1.10 

708 

275 

it  is  almost  certain  that  a  combination  of  the  steeply-sloped  and  more 

shallowly-sloped  spectral  components  was  being  measured,  but  because  of 

Che  shorter  detrend  interval,  we  were  still  primarily  sensitive  to 

the  shallowly-sloped  component-  Thus,  our  distribution  for  p  Is 

B 

shallower  than  would  be  expected  from  the  Index  values  we  reported  In 
our  earlier  summaries  F Rlno.  1979],  although,  as  we  have  already  noted, 
noise  contamination  Is  introducing  small  biases  in  the  results.  Also, 
Fremouw  et  al.  [1985a]  have  noted  that  in  the  early  analysis  of  the 
Wideband  data  detrending  without  subsequent  windowing  was  used.  This 
allowed  enough  sldelobe  leakage  to  bias  the  results  toward  more  shallow 
values . 

To  substantiate  the  general  features  of  the  spectral  character¬ 
istics,  we  note  that  processes  that  are  unique  to  the  auroral  zone 
should  show  a  dependence  on  invariant  latitude,  most  likely  controlled 
by  magnetic  activity.  Thus,  In  Figure  15  we  have  plotted  the  single¬ 
component  p  Index  values  and  the  p.  and  p_  Indices  for  the  dominant  Type 

A  B 

3  spectra.  It  can  be  seen  that  both  the  p  and  the  p.  Indices  tend  to 

A 

decrease  with  Increasing  Invariant  latitude;  moreover,  because  the  trend 

is  asymmetric.  It  is  most  likely  to  be  a  true  latitudinal  dependence 

rather  than  an  elevation-angle  or  magnetic  aspect-angle  dependence.  The 

corresponding  data  for  p  ,  by  comparison,  show  a  distinct  tendency  to 

B 

form  a  local  maximum  at  the  point  where  the  propagation  path  Intercepts 
an  L  shell.  This  effect  has  been  discussed  by  Fremouw  et  al .  [1985a]. 
who  attribute  It  to  a  scale-size  dependent  anisotropy.  The  low- 
frequency  enhancement,  however,  appears  to  be  a  genuine  auroral-zone 
characteristic. 


15(a).  Invariant  latitude  dependence  of  auroral-zone  spectral  indices --Poker  Flat, 
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DISCUSSION 

The  data  analysis  presented  In  this  paper  provides  a  consistent 
average  spectral  characterization  of  F-region  Irregularities  having 
scale  sizes  from  a  few  tens  of  k.ilometers  to  a  few  hundreds  of  meters. 
Within  this  scale-size  regime,  equatorial  irregularities  admit  one  of 
two  characterizations:  (1)  single  power  laws  having  one-dimensional  in 
situ  spectral  indices  near  but  somewhat  larger  than  2.  and  (2)  two- 
component  power  laws  having  intermediate  scale  indices  less  than  2  with 
a  transition  between  500  and  800  m  to  power  laws  more  steeply  sloped 
than  2.  From  the  scintillation  data  alone,  we  are  not  able  to  identify 
the  conditions  under  which  each  type  occurs;  however,  in  situ  rocket 
probes  have  shown  that  the  single-component  power  laws  are  observed 
primarily  in  bottom-side  structures  [Keskinen  et  al..  1981]  and  within 
or  on  the  walls  of  large  depletions  F Rlno  et  al . .  1981a;  Kellev.  et  al . . 
1982;  Kellev  et  al . .  1986;  LaBelle  and  Kellev.  1986].  In  situ  satellite 
probes,  which  scan  the  structures  horizontally  rather  than  vertically, 
seem  to  be  in  agreement  with  this  occurrence  pattern  of  single-component 
and  two-component  spectra  F Livingston,  private  communication. 

The  single-component  spectra  having  one-dimensional  in  situ  indices 
near  2  agree  with  the  predictions  of  both  theory  and  numerical  simula¬ 
tions  of  the  gravity-driven  Rayleigh-Taylor  instability  F Ossakow.  1979; 
Keskinen  et  al. .  1981];  however,  as  noted  by  Kellev  et  al.  [1986],  the 
two-component  spectra  remain  largely  unexplained.  LaBelle  et  al.  [1986] 
have  reviewed  two  hypptheses  that  can  possibly  explain  the  two-component 
spectra.  One  hypothesis  assumes  that  the  structure  is  caused  by  two- 
dimensional  turbulence  in  the  velocity  field,  which  can  give  rise  to  a 
two-component  spectrum  in  connection  with  a  reverse  cascade  of  energy 
(from  small-  to  large-scale  sizes),  and  a  downward  cascade  of  enstropy 


originally  proposed  by  Kraichnan  [1967];  it  has  been  observed  in  neutral 
turbulence  F Lilly.  1983]  and  has  been  demonstrated  in  numerical  simu¬ 
lations  of  "guiding-center  plasmas"  by  Sevier  et  al.  [1975].  As 
reviewed  in  Klntner  and  Sevier  [1985],  however,  the  application  of  the 
theory  relies  on  a  passive  scalar  relation  between  density  and  velocity 
(collisionless  structuring)  and  it  implicitly  requires  an  isotropic 
velocity  field.  We  believe  that  these  restrictions  are  too  severe  and 
that  alternate  hypotheses  must  be  considered. 

We  had  suggested  an  alternative  hypothesis,  which  is  discussed  in 
LaBelle  et  al.  [1986].  Batchelor  [1959]  showed  that  highly  anisotropic 
structures  in  a  passive  scalar  field  give  rise  to  a  h,  ^  spectrum. 
Batchelor  was  studying  passive  scalar  turbulence  beyond  the  viscous  cut¬ 
off  in  the  velocity  field  where  strong  linear  shears  are  known  to  de¬ 
velop.  We  noted  that  a  very  similar  structure  occurs  in  the  velocity 
field  of  the  highly  anisotropic  evolution  of  convective  instabilities. 
The  irregularities  are  drawn  out  along  the  direction  of  the  average 
cross-field  drift,  which  may  be  caused  by  any  combination  of  electric 
fields,  gravity,  or  neutral  winds.  The  results  show  that  a  very 
shallowly-sloped  intermediate-scale  regime  is  Indeed  associated  with  the 
highly  anisotropic  structure  evolution  in  purely  collislonal  instabili¬ 
ties;  however,  the  results  also  suggest  that  the  anisotropy  is  scale- 
size  dependent.  In  any  case,  the  two-component  model  itself  seems  to  be 
well  supported  by  a  variety  of  Independent  data — if  not  yet  fully 
explained . 

Regarding  the  steeply-sloped  transition-scale  regime,  the  spectral 
Indices  derived  from  the  phase  scintillation  data  give  results  signif¬ 
icantly  shallower  than  values  reported  from  in  situ  measurements  and 
even  intensity  scintillation  data  FKellev  et  al. .  1982;  Franke  and  Liu. 
1983].  As  discussed  in  Report  I,  it  is  possible  that  the  diffraction- 
induced  rapid  phase  changes  are  yet  biasing  the  data.  Indeed,  diffrac¬ 
tion  effects  in  phase  tend  to  become  larger  as  the  scale  sizes  decrease, 
whereas  intensity  spectra  are  affected  very  differently  and  would  tend 


to  give  a  steeper  result  than  is  measured  in  situ.  Also,  at  suffi¬ 
ciently  small  scale  sizes,  temporal  rearrangement  of  the  structure 
during  the  measurement  can  also  bias  or  smear  the  result.  In  any  case, 
it  is  clear  that  the  steepened  portion  of  the  spectrum  is  a  manifesta¬ 
tion  of  enhanced  diffusion-like  processes  that  are  rapidly  eroding  the 
steep  gradients  that  would  otherwise  result  from  the  natural  evolution 
of  convective  structuring. 

Turning  to  the  auroral-zone  data,  the  existence  of  a  low-frequency 
spectral  enhancement  in  t.ie  Wideband  satellite  data  had  been  noted  by 
Fremouw  et  al..  [1985a];  however,  the  usual  practice  of  fitting  linear 
segments  over  fixed  temporal-frequency  regimes  seems  to  have  obscured 
the  persistence  of  this  feature.  Indeed,  as  with  the  equatorial  Type  2 
spectra,  the  occurrence  of  auroral-zone  Type  3  spectra  is  remarkably 
consistent  over  time  and  auroral-zone  location.  Auroral-zone  irregu¬ 
larities  have  been  reviewed  by  Feier  and  Kellev  [1980],  Kellev  et  al . . 
[1982],  Vickrey  and  Kellev  [1982],  and  Tsunoda  [1986],  who  discussed  a 
general  framework  for  the  production,  convection,  and  dissipation  of 
auroral-zone  irregularities.  They  emphasized  the  connection  of  the 
hlgh-latltude  F  region  to  the  distant  magnetosphere  and  to  a  highly 
conducting  E  region,  which  can  provide  .direct  sources  of  irregularities 
through  structured  velocity  fields  and  precipitating  particles,  and  can 
profoundly  influence  the  evolution  of  such  structures. 

It  is  possible  that  the  low-frequency  component  of  the  two- 
component  spectra  is  a  manifestation  of  direct  velocity  stirring — as 
suggested  by  Vickrev  and  Kellev  [1982].  Once  established,  convective 
instabilities  would  sustain  a  continued  structure  enhancement.  In  this 
context,  the  patchiness  of  high-latitude  scintillation  would  be  attrib¬ 
uted  to  electrodynamic  processes  that  cause  structuring  on  a  scale  of 
tens  to  hundreds  of  kilometers  f Robinson  et  al . .  1985;  Tsunoda  et  al.. 
1985;  Tsunoda.  1986].  Fremouw  et  al..  [1985a]  have  suggested  that  the 
sheet-like  anisotropy  observed  in  the  night-side  auroral-zone  scintil¬ 
lation  data  is  also  caused  by  the  direct  deposition  of  structure.  With 
strong  zonal  convection,  the  structures  would  presumably  be  Imparted 


45 


with  a  configuration  close  to  the  form  toward  which  any  east-west 
gradient  would  tend  to  evolve,  as  discussed  by  Vlclcrev  and  Kelly 
[1982]. 


The  source  of  the  more  shallowly-sloped  component  beyond  one 
kilometer  is  puzzling.  The  scale  sizes  are  small  enough  that  diffusion 
should  be  acting  to  steepen  the  spectrum.  We  note  that  Villain  et  al. 
[1986]  report  a  steepening  of  in  situ  spectra  beyond  a  few  hundred 
meters.  This  suggests  that  the  shallowly-sloped  portion  at  the  high- 
frequency  end  of  the  spectrum  might  well  be  caused  by  the  same  mechanism 
that  causes  shallowly-sloped  spectra  In  the  equatorial  zone,  but  extend¬ 
ing  to  somewhat  smaller  scale  sizes.  The  vastly  different  geometries 
causing  variations  along  (rather  than  across)  field  lines  could  possibly 
account  for  such  effects.  Finally,  we  note  that  Basu  et  al.  [1984]  have 
observed  shallow-sloped  in  situ  spectra  to  occur  in  the  vicinity  of 
strong  velocity  shears  and  have  suggested  a  cause-effect  relation.  All 
the  analyses  performed  to  date,  however,  show  that  inertial  effects 
alone  are  either  stabilizing  or  that  they  evolve  structure  at  a  much 


slower  rate  than  purely  colllslonal  mechanisms 


[.  1984; 


..  1985; 


..  1986]. 


Ultimately,  a  consistent  explanation  of  the  results  will  require  a 
better  understanding  of  how  intermediate-scale  structures  evolve,  and 
particularly  how  this  evolution  manifests  itself  in  the  measurable  one- 
dimensional  spectral  characteristics.  Tsunoda  [1986]  has  developed  a 
unified  framework  within  which  the  auroral-zone  structure  enhancement 
can  be  Interpreted.  However,  much  of  the  observed  differences  between 
the  equatorial  and  auroral-zone  structures,  may  yet  be  explainable  by 
the  different  geometries  or,  simply,  by  the  different  conductivity 
ratios . 

As  was  described  in  Rep»ort  I,  the  propagation  results  reported  here 
must  ultimately  be  extrapolated  to  predict  the  nuclear  disturbed  environ¬ 
ment.  The  first  link  in  that  extrapolation  procedure  Involves  verifying 
that  the  phase  spectral  shapes  reflect  the  in  situ  density  structure. 


This  can  only  be  done  by  conducting  an  analysis  similar  to  this  one  on 
in  situ  measurements.  Such  an  effort  is  currently  underway  for 
Atmosphere  Explorer  equatorial  data.  We  intend,  under  continued  DNA 
support,  to  analyze  high-latitude  in  situ  data  as  well.  In  addition  to 
verifying  the  spectral  shapes,  these  studies  will  emphasize  the  physics 
that  control  the  different  shapes  observed.  It  is  only  through  a  de¬ 
tailed  understanding  of  these  underlying  physical  processes  that  the 
extrapolation  to  the  nuclear  regime  can  be  done  with  confidence. 
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APPENDIX 


AN  ALGORITHM  FOR  SEGMENTED  LINEAR-CURVE  FITS* 


A.l  INTRODUCTION. 

We  have  developed  an  automated  procedure  for  malting  multisegmented 
fits  to  functional  data.  In  the  past,  it  was  necessary  to  fit  linear 
segments  to  data  by  malting  linear  least-squares  fits  to  specific  regions 
of  the  data.  The  ranges  for  the  linear  fits  are  preselected  or  the 
linear  regions  in  the  data  are  determined  by  observing  each  data  set 
individually.  While  the  former  method  is  very  efficient,  it  is  sub¬ 
jective;  furthermore,  it  cannot  be  used  if  the  linear  regions  are  non- 
statlonary.  The  latter  method,  while  objective,  is  very  time  consuming, 
thus.  It  cannot  be  used  to  analyze  a  large  data  set.  The  new  algorithm 
combines  objectivity  with  efficiency;  it  is  a  superior  method  for  sort¬ 
ing  data  into  linear  regimes.  In  this  Appendix,  we  describe  the  method 
and  its  usage;  then  we  explain  the  details  of  the  algorithm. 


A. 2  PHILOSOPHY. 

Like  the  method  of  least  squares,  our  scheme  approximates  the 
dependence  of  one  variable  on  another.  Both  methods  describe  the 
relationship  analytically.  The  method  of  least  squares  finds  a 
polynomial  representation  of  the  form 


F(x)  =  Ajj  •  x"  ^  f  •  x"  ^  +  ...  +  ^2  * 


This  Appendix  has  been  prepared  by  Terl  M.  Dabbs 
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^  ^  *  m  *  _M  *  ft  ^  ft  ^  M  ^  ■  * 


Our  scheme  finds  a  mult isegmented  representation  of  the  form 


F(X)  =  •  X  +  .  X^  <  X  <  X^. 

=  «2  ’  ^  ®2  .  ^  X  <  X3. 


X  +  B  ,  .  X  <  X  <  X 
n-1  n-l  n 


Although  the  two  methods  describe  the  relationship  between  two  variables 
in  a  similar  fashion,  each  representation  has  its  own  advantages  and 
disadvantages.  We  developed  the  multlsegmented  fitting  scheme  to  com¬ 
pensate  for  the  shortcomings  in  the  method  of  least  squares.  To  be 
specific,  we  wanted  a  scheme  that  would  fit  two-dimensional  data  with 
the  flexibility  of  a  hlgh-order-polynomlal  fit,  and  that  would  give 
results  that  are  as  simple  to  Interpret  as  a  linear  fit. 

Every  two-dimensional  data  set  can  be  accurately  represented  by 
some  high-order  polynomial;  when  the  degree  of  the  polynomial  fit  is  one 
less  than  the  number  of  data  points,  the  pxjlynomlal  can  represent  the 
data  exactly.  Unfortunately,  the  significant  structure  in  the  data  is 
not  directly  Inferred  from  the  high-order-polynomial  fits.  Although  the 
continuous  changes  in  the  polynomial  can  be  found  Indirectly  by  differ¬ 
entiating  the  equation,  the  local  trends  are  still  unavailable;  thus, 
comparisons  between  polynomials  are  complicated  and  difficult  to  com¬ 
prehend.  In  contrast,  the  trend  in  a  segment  is  obvious  from  its  linear 
equation.  The  slope  of  the  line  defines  the  trend  in  the  data  over  the 
range  of  the  segment.  Because  the  slof)e  of  the  segment  measures  the 
trend  in  the  data,  comparisons  between  segments  are  direct  and  simple  to 
understand. 

Although  linear  fits  are  easy  to  Interpret,  they  will  not  accu¬ 
rately  represent  trends  in  complicated  data.  They  will  hide  significant 
structure  in  data  when  their  range  covers  two  or  more  distinct  regions. 
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For  example,  the  data  in  Figure  16  have  two  distinct  regions.  For  X  <  0, 
the  values  of  Y  are  increasing;  for  X  >  0,  the  values  of  Y  are  decreas¬ 
ing.  A  linear  least-squares  fit  approximates  the  data  by  the  linear 
equation  F(X)  =  0.5  *  X  +  1.0.  The  slope  of  the  line,  0.5,  shows  a 
positive  dependence  of  Y  on  X.  The  two  significant  trends  have  been 
averaged.  If  the  ranges  of  the  two  trends  could  be  defined,  the  trends 
could  be  separated  by  malting  a  linear  fit  to  each  regime.  For  Instance, 
if  we  made  two  linear  fits  to  the  data  in  Figure  16,  a  fit  less  than 
zero  and  another  fit  greater  than  zero,  we  would  measure  the  two  trends. 
The  first  region  would  be  approximated  by  the  equation 

F(X)  =  2.0  •  X  +  5.0 

and  the  second  region  would  be  approximated  by  the  equation 

F(X)  =  -1.0  •  X  +  5.0 

The  two  slopes,  2  and  -1,  accurately  describe  the  relationship  between  X 
and  Y  for  the  two  regions. 

Clearly,  when  we  Itnow  the  ranges  of  the  linear  regions  in  the 
data,  we  can  give  a  simple  and  accurate  description  by  malting  linear 
least-squares  fits  to  each  linear  regime.  Unfortunately,  we  do  not 
always  Itnow  their  ranges.  The  linear  regions  can  be  randomly  varying — 
they  can  move  around  in  the  data.  For  example,  the  spectra  of  phase 
scintillation  data  admit  both  random  and  systematic  variations.  A 
beacon  on  a  satellite  transmits  radio  frequencies  that  measure  phase 
scintillation.  Because  the  geometry  of  the  satellite  pass  changes  the 
effective  speed  of  the  satellite  through  the  ionosphere,  the  relation¬ 
ship  between  the  spatial  domain  and  the  frequency  domain  changes. 
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Assuming  the  linear  regions  in  the  phase  scintillation  data  are  defined 
spatially,  wo  know  their  positions  in  frequency  will  vary;  therefore,  we 
cannot  define  the  ranges  of  the  linear  regions  in  terms  of  the 
frequency. 

Because  we  cannot  always  preselect  the  linear  regions  in  func¬ 
tional  data,  we  developed  a  scheme  that  uses  the  data  to  find  their  own 
linear  regimes.  By  predicting  the  ranges  of  the  linear  regions,  we  are 
able  to  make  multlsegmented  fits.  These  fits  are  simple  to  Interpret 
and  they  give  an  accurate  description  of  the  data.  This  method  can  be 
used  whenever  a  least-squares  polynomial  fit  would  be  an  appropriate 
approximation  to  the  data.  We  use  this  scheme  to  sort  the  data  into 
linear  regimes,  when  the  ranges  of  the  linear  regions  cannot  be  pre¬ 
determined. 

A. 3  DESCRIPTION  OF  MULT I SEGMENTED  FITTING  PROCEDURE. 

By  observing  that  different  people  tend  to  find  the  same  linear 
segments  in  data  when  they  are  presented  in  a  rectangular  graph,  we 
concluded  that  the  linear  regimes  are  selected  objectively  and  system¬ 
atically.  Furthermore,  we  believed  that  an  algorithm  could  emulate  this 
process.  By  decomposing  this  abstract  visual  experience  into  simple 
tasks,  we  develop>ed  a  logical  procedure  for  finding  linear  trends  in 
two-dimensional  data.  The  resulting  algorithm  consists  of  finding  the 
overall  trend,  determining  the  linear  regions,  then  finding  the  linear 
equations  that  represent  each  region. 

We  find  the  overall  trend  by  making  a  least-squares  polynomial  fit 
to  the  data,  then  evaluating  the  polynomial  at  the  end  points.  The  line 
containing  these  two  points  mathematically  defines  the  overall  trend — it 
describes  the  change  in  the  dependent  variable  over  the  range  of  the  in¬ 
dependent  variable.  We  have  not  attempted  to  make  this  scheme  general 
enough  to  determine  the  ideal  degree  of  the  polynomial  fit;  instead,  we 
let  the  person  using  this  scheme  choose  the  degree  of  the  fit.  The 
polynomial  should  be  flexible  enough  to  bend  with  the  significant  trends 
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in  the  data;  thus,  the  degree  of  the  polynomial  should  be  greater  than 
or  equal  to  the  number  of  linear  regimes.  Figure  17  shows  this  trend. 


From  the  polynomial  fit  and  the  overall  trend,  we  estimate  the 
linear  regions,  i.e.,  the  regions  of  the  polynomial  that  are  monotoni- 
cally  increasing  or  decreasing  relative  to  the  overall  trend.  The  end 
points  of  each  linear  region  have  a  unique  relationship  to  the  overall 
trend.  At  these  points,  the  tangent  to  the  polynomial  is  parallel  to 
the  overall  trend.  In  mathematical  terms,  let 

F  =  least-squares  polynomial  fit  to  the  data 

F'  =  derivative  of  least-squares  polynomial  fit  to  the  data 


1  ■, 

f: 


Slope  =  slope  of  line  defining  overall  trend. 

We  find  the  points  within  the  range  of  the  data  where 


F'  =  slope 

That  is,  we  find  the  solutions  to  the  following  equation. 


I  F'  -  slope  =  0 


The  degree  of  the  equation  determines  the  maximum  number  of  solutions  to 
this  equation.  Because  we  are  using  these  solutions  to  separate  the 
data  Into  increasing  and  decreasing  regions,  we  are  only  interested  in 
the  solutions  that  occur  between  the  end  points.  The  mean  value  theorem 
ensures  that  at  least  one  solution  will  occur  in  this  range.  Using  the 
end  points  and  the  solutions  between  them,  we  define  the  regions  in  the 
polynomial  that  are  monotonlc  relative  to  the  overall  trend.  Because 
data  are  discrete  rather  than  continuous,  we  must  check,  that  these 
monotonlc  regions  are  definable.  When  the  degree  of  a  polynomial  fit  is 
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higher  than  the  data  require,  the  polynomial  can  make  bends  between  two 
consecutive  data  points.  These  bends  are  meaningless;  they  do  not 
represent  the  data.  Again,  we  have  not  attempted  to  make  this  scheme 
general  enough  to  determine  the  number  of  data  points  that  would  make  a 
region  significant.  Because  a  segment  cannot  be  defined  without  two 
distinct  points,  each  region  must  contain  at  least  two  data  points; 
however,  the  person  who  uses  this  scheme  may  define  a  more  stringent 
requirement . 

If  any  region  does  not  meet  this  criterion,  we  make  a  new  poly¬ 
nomial  fit  to  the  data.  The  degree  of  the  new  polynomial  fit  Is  one 
less  than  the  previous  fit;  thus,  the  new  fit  Is  less  flexible.  We 
repeat  the  entire  process  until  we  find  monotonic  regions  in  the 
polynomial  that  are  representative  of  the  data.  Figure  18  shows  the 
points  that  estimate  the  linear  regions  in  the  data.  Notice  that  the 
tangents  to  the  polynomial  at  these  points  are  parallel  to  the  overall 
trend. 

After  predicting  the  linear  ranges  in  the  data,  we  refine  our 
estimates  by  making  linear  least-squares  fits  to  each  regime.  If  the 
data  are  well-behaved,  we  make  the  linear  fits  to  the  original  data 
points;  otherwise,  we  make  the  fits  to  the  polynomial  representation  of 
the  data.  Figure  19  illustrates  the  linear  fits  to  the  original  data 
points  over  each  region,  and  Figure  20  Illustrates  the  linear  fits  to 
the  polynomial  representations.  Although  the  fits  to  the  original  data 
points  are  more  representative  of  the  data,  these  fits  can  be  adversely 
affected  by  stray  data  points.  From  these  linear  fits  we  redefine  the 
ranges  of  each  linear  region.  We  determine  the  corrected  range  of  each 
region  by  finding  the  intersection  of  the  lines  representing  consecutive 
regions.  Figures  21  and  22  show  the  final  segment  fits  to  the  data,  and 
to  the  polynomial  representation  of  the  data,  respectively.  Although 
Figure  21  is  a  better  fit,  both  fits  emphasize  the  major  trends. 


data  over  each  estimated  linear  region. 


A. 4  CONCLUSION. 


We  have  developed  a  unique  scheme  for  sorting  two-dimensional  data 
into  linear  regimes.  Unlike  other  methods  that  fit  segments  to  data, 
this  algorithm  determines  the  number  of  linear  regimes  in  the  data  and 
the  ranges  of  these  regions.  Because  it  determines  the  locations  of 
each  region  from  a  polynomial  model,  the  accuracy  of  the  linear  fits  is 
dependent  upon  the  polynomial's  ability  to  represent  the  data.  When  the 
polynomial  fit  accurately  represents  the  data,  the  linear  fits  will  also 
accurately  represent  them.  This  scheme  provides  an  efficient  and 
objective  way  to  make  multisegmented  fits  to  two-dimensional  data. 
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